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AERONAUTIC SYMBOLS 


1. FUNDAMENTAL AND DERIVED UNITS 



Symbol 

Metric 

English 

Unit 

Abbrevia- 

tion^ 

i 

Unit 

Abbrevia- 

tion 

Length _ _ _ _ 

l 

meter 

m 

foot (or mile) 

ft (or mi) 

Time 

t 

second 

s 

second (or hour) 

sec (or hr) 

Force 

F 

weight of 1 kilogram * 

kg 

weight of 1 pound 

lb 

Power 

P 

horsepower (metric) 


horsepower 

hp 


TT 

/kilometers per hour 

kph 

miles per hour _ _ 

mph 


V 

\ meters per second. _ _ 

mps 

feet per second. 

fps 


2. GENERAL SYMBOLS 


Weight = mg 

Standard acceleration of gravity =9.80665 m/s 3 
or 32.1740 ft/sec 3 

Mass=— 

9 

Moment of inertia= «■&&*. (Indicate axis of 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


v Kinematic viscosity 

p Density (mass per unit volume) 

Standard density of dry air^ 0.12497 kg-m“ 4 -s 9 at 15° O 
and 760 mm; or 0.002378 lb-ft“ 4 sec 2 
Specific weight of “stan dard ”?, air, 1.2255 kg/m 3 or 
0.07651 lb/cu ft 


3. AERODYNAMIC SYMBOLS 


Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio, ^ 

True air speed 
Dynamic pressure, ^pF 3 

Lift, absolute coefficient 0 L = ^ 

Drag, absolute coefficient 

Profile drag, absolute coefficient 

Induced drag, absolute coefficient C D .=Q 

Parasite drag, absolute coefficient C °>=fs 

Cross-wind force, absolute coefficient C c —z 


_A 


lio 

it 

Q 

R 


a 

€ 

«0 

OLi 

Ota 


Angle of setting of wings (relative to thrust line) 
Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 

Reynolds number, p— where l is a linear dimen- 

sion (e.g., for an airfoil of 1 .0 ft chord, 100 mph, 
standard pressure at 15° C, the corresponding 
Reynolds number is 935,400; or for an airfoil 
of 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 

Flight-path angle 
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THE USE OF SOURCE-SINK AND DOUBLET DISTRIBUTIONS EXTENDED TO THE SOLUTION 
OF BOUNDARY-VALUE PROBLEMS IN SUPERSONIC FLOW 

By Max. A. Heaslet and Harvard Lomax 


SUMMARY 

A direct analogy is established between the use of source-sink 
and doublet distributions in the solution of specific boundary- 
value problems in subsonic wing theory and the corresponding 
problems in supersonic theory . The concept of the 11 finite part " 
of an integral is introduced and used in the calculation of the 
improper integrals associated with supersonic doublet distri- 
butions . The general equations developed are shoum to include 
several previously published results and particular examples 
are given for the loading on rolling and pitching triangular 
wings with supersonic leading edges. 

INTRODUCTION 

The problem of finding pressure distributions over airfoils 
of arbitrary shape and plan form or of finding airfoils which 
have arbitrary pressure distributions is one of the most 
fundamental problems in aerodynamic theory. At the pres- 
ent time the most important and satisfactory approach to 
problems of this type is provided by the methods of so-called 
thin-airfoil theory. The essential assumptions in this theory 
are that the perturbation velocities induced by the airfoil 
are small relative to the free-stream velocity and that the 
boundary conditions can be specified in a fixed reference 
plane. 

Under the assumptions of thin-airfoil theory the theoreti- 
cal analysis of a problem in wing theory resolves itself into 
the task of determining the solution of a second-order linear 
partial differential equation with prescribed boundary 
conditions. In the case of purely subsonic flow, Laplace's 
equation in three dimensions must be considered, while in 
purely supersonic flow the differential equation which arises 
is algebraically equivalent to the two-dimensional wave 
equation of mathematical physics. The classical solutions 
of these two equations have been developed along two dis- 
tinct lines: First, by use of orthogonal functions which can 
be derived in terms of the boundary conditions, and alter- 
natively by means of Green's theorem which in turn utilizes 
a known particular solution of the partial differential equa- 
tion together with the given boundary conditions. 

One particular solution associated with Laplace's equation 
and subsonic aerodynamics has been found to be outstanding 
in its mathematical usefulness and, when identified with the 
velocity potential, has a physical interpretation which can 
supply, in direct application, added insight into the nature of 
the problem. This function is referred to as the “funda- 


mental solution" and can be developed from the concept of 
a so-called source. A concomitant development to the 
source potential is the doublet potential, and appropriate 
distributions of these functions are known to be sufficient 
for the solution of all problems in subsonic wing theory. 

The extension of the use of the fundamental solution to 
problems in purely supersonic flow introduces mathematical 
difficulties which differ essentially from those encountered 
at low speeds. Both the source and the doublet potentials 
possess singularities on their conical characteristic surfaces 
or Mach cones and, in the case of the doublet, the singularity 
is of higher order than can be treated by elementary mathe- 
matical methods. In the historical development of the 
solutions of the wave equation this trouble was circumvented 
by replacing the source potential by other particular solu- 
tions of the differential equation. As an example, Volterra 
(reference 1) introduced the integral of the fundamental 
solution and in that way reduced the order of the singularities 
involved. The analytical development of Volterra's theory 
presents no inherent difficulties (e. g., reference 2) but the 
physical significance of the particular solution is lost, the 
direct analogy with subsonic theory no longer exists, and a 
certain amount of mathematical inefficiency arises since, 
after using the integral of the source potential, it is found 
necessary to resort at the end of the analysis to taking a 
final derivative. 

In this report, following methods introduced by Hadamard 
(reference 3), a general solution to the thin-airfoil problem 
in supersonic theory will be given in terms of the distribution 
of sources and doublets over the given reference plane. 
Furthermore, a discussion of the nature of the boundary 
values required will be given. For properly set problems 
in wave theory it has been found necessary to specify, 
usually, both the required function and its derivative with 
respect to time along the boundary considered. In aero- 
dynamic applications of the wave equation associated with 
lifting-surface theory and thickness distributions it will be 
shown that only a knowledge of the unknown function or 
its normal derivative along the boundary is needed since a 
relationship between the two functions will be established 
on the boundary surface. 

In the theoretical portion of the report a brief presenta- 
tion will be made of the differential equations involved and 
the two forms of the fundamental solution. An outline is 
then given of the types of boundary-value problems encoun- 
tered and, since the purpose of the report is to extend the 
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concepts of thin-airfoil theory which are used in subsonic 
theory to problems arising in supersonic theory, a discussion 
will be given of the subsonic development as a basis for the 
analogy which exists between the methods of solution corre- 
sponding to the two regimes of flow. In the discussion of 
the purely supersonic case it will be shown that the intro- 
duction of the concept of “finite part” will provide a tech- 
nique whereby the improper integrals arising from the use 
of doublets may be evaluated in a straight-forward manner. 
The applications of the theoretical developments will include 
the rederivation of some previously published results and 
will also contain the calculation of load distributions for 
rolling and pitching triangular wings with leading edges 
swept ahead of the Mach cone from the vertex of the 
triangles. 

SYMBOLS 


b 

c 

M 

n 

n u n 2) n z 

V 

P 

2 

Q 

r 

T s 
To 

R 

S 


V 

X,Y,Z 
x, y, z 
|3 

e 

X 

v 

Vl, ”2, U3 
a 


T 

$ 

ft 


c v 


span of wing 
chord of wing 
free-stream Mach number 
normal to arbitrary surface 
direction cosines of normal n 
static pressure 
rate of roll about A" axis 
free-stream dynamic pressure 
rate of pitch about Y axis 
V(x— xQ 2 +(y— yi) 2 +(z— zQ 2 
V(x z -x 1 ) 2 +(y— y i) 2 +g~ 2 
V (s— xQ 2 — (y— yQ 2 — (z — zQ 2 
V(x— x : ) 2 — (y— y,) 2 — s 2 
arbitrary region of integration 
surface enclosing region R 

perturbation velocities in direction of X , Y, and Z 
axes respectively 
free-stream velocity 

Cartesian coordinates in original space variables 
transformed system of coordinates 

VIM^I 

infinitesimal used in analysis 

surface along which stream enters induced field of 
wing 

conormal to arbitrary surface 
direction cosines of conormal 

variable representing either acceleration potential, 
velocity potential, or any of the three perturba- 
tion velocity components 
surface on which boundary conditions are given 
perturbation velocity potential 
variable representing either acceleration potential, 
velocity potential, or any of the three perturba- 
tion velocity components 
pressure coefficient 


Aj) 

T 


load coefficient 


O x 

c. 


V 2 

□ 2 

r 

u 


differential operator (^2+^2+^) 

, /d 2 d 2 d 2 \ 

differential operator ) 

sign denoting “finite part” of integral 


SUBSCRIPTS 

subscript denoting value of variable on upper sur- 
face of wing 

l subscript denoting value of variable on lower sur- 

face of wing 

j subscript denoting variable of integration 

c subscript on r denoting fundamental solution in 

supersonic flow 

SUPERSCRIPT 

' superscript denoting value of variable on opposite 

side of t from fixed point (x, y , z) 

THEORETICAL DEVELOPMENT 

LINEARIZED EQUATIONS AND BOUNDARY CONDITIONS 

The linearization of the second-order differential equation 
for compressible fluid flow is developed under the assump- 
tions of thin-airfoil or small-perturbation theory. If the 
velocity vector of the free stream is parallel to and in the 
direction of the positive X axis, the resulting differential 
equation is expressible in the form 


„ Tijf2\ d 2 0 ,c> 2 2 c«2_ n 

i\-M ) bX 2+ 572+52 2-0 


( 1 ) 


where ft represents a^velocity potential, acceleration poten- 
tial, or any one of the perturbation velocities while M is the 
constant value of the free-stream Mach number. Assuming 
the plane of symmetry of the airfoil to lie in the X Y plane, 
the boundary conditions associated with equation (1) are 
given for Z= 0. Moreover, if u,t;,and w represent, respective- 
ly, the perturbation velocity components along the X, Y, 
and 7 axes, and if the velocity of the free stream is V, the 
direction cosines of any streamline are proportional to the 
point functions V-\-u(X, Y, Z), v(X, Y, Z), and w(X } Y, Z) 
while pressure coefficient C v is given by the relation 


r -_?2f 

Op— y 


( 2 ) 


... rr . /moment about X axis\ 

rolling-moment coefficient ( , • ) 

V go X wing area / 

d<7i 

d(Pb/2V ) 


Detailed discussions of these results may be found in refer- 
ence 4. 

Introducing the affine transformations 

y=V±( i-M 2 )F 


z= V± (1— M 2 )Z 

where the signs under the radicals are chosen so that real 
values result, it follows that in the subsonic case (M< 1) 
equation (1) reduces to 


^ft 

dx 2 ^~dy 2 ^dz 2 U 


( 3 ) 
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while the supersonic case (M^> 1) yields 


d 2 ^_d 2 ^l_d 2 9. = 
dx 2 dy 2 dz 2 

The fundamental solution associated with equation (3) 


(4) 

is 


I=[(z-Zi ( 2 +(2/-2 /i) 2 +0-Zi) 2 ] 2 (5) 

or, in terms of the original space variables, 

I=[(X-X0 2 +/3 2 (y-r0 2 +/J 2 (Z-Z0 2 U (5a) 


where 

/3 2 =(1-M 2 ) 

When the wave equation is to be considered the fundamental 
solution takes the form 


y = [(x-x 1 ) 2 -(y-yi) 2 -(z-z l ) 2 ] 2 (6) 

or 

} 1 =[(X-X 1 ) 2 -^ 2 (F-F 1 ) 2 -^ 2 (Z-Z 1 ) 2 ]" i (6a) 

where 

/3 2 = (M 2 — 1) 

These fundamental solutions are directly related both in 
subsonic and supersonic flow to the velocity potentials at 
the point (x, y, z) or ( X , Y, Z ) of unit sources situated at the 
point (a?i, yi, z x ) or (X ly Ti, Z x ). The velocity potential of a 
doublet may be obtained by taking a directional derivative 
of the source potential, the direction of the axis of the doublet 
coinciding with the direction along which the derivative is 
taken. These two functions will be seen to be of paramount 
importance when Green’s theorem is applied to the given 
boundary conditions. 

It remains now to mention the types of boundary condi- 
tions which appear in problems associated with wing theory. 
As a convenience to the development of the theory the 
normalized forms (equations (3) and (4)) of equation (1) 
will be used and boundary conditions will be assumed known 
with respect to the x , y , 2 coordinate system. Retrans- 
formation to the X , F, Z system of axes can be made quite 
simply wherever needed in application. In order to define 
the boundary conditions, two subscripts will be introduced: 
The first, u, denotes the value of the required function on the 
upper surface, that is, the limit of the function as 2 ap- 
proaches zero from the positive direction; the second, Z, 
denotes the value on the lower surface, that is, the limit of the 
required function as 2 approaches zero from the negative 
direction. 

Using these definitions the four boundary-value problems 
of principal interest can be defined as follows: 

1 . Symmetrical nonlifting airfoil with specified slope. — 
In this case w u =w t =0 over all of the xy plane except for 
region occupied by the airfoil where 2w u = — 2w t = Aw=f(x, y), 
the function being determined by the geometry of the wing. 
Over all of the xy plane, Au=0. 


2. Lifting plate with specified loading.— It is given that 
Au=u u —Ui=0 over the xy plane except for the region 
occupied by the airfoil where Au=f(x, y), the function being 
determined by the specified loading. Moreover, Aw=0 
everywhere. 

3. Lifting plate with specified camber, twist, and angle of 
incidence. — Over the xy plane Aw= 0 everywhere. And, 
except for the region occupied by the airfoil, A / m= 0. Over 
the region occupied by the airfoil w=j{x , y) where /(a;, y) is 
determined by the given camber, twist, and angle of 
incidence. 

4. Symmetrical nonlifting airfoil with specified pressure 
distribution. — Over the xy plane Au=0 everywhere. And, 
except for the region occupied by the airfoil Aw=0. Over 
the region occupied by the airfoil C p =f(x, y) where /(x, y) is 
given. 

It should be pointed out that the first two problems con- 
sidered here differ from the usual type of boundary-value 
problem encountered. In the so-called Dirichlet or Neu- 
mann problems, which arise in connection with Laplace’s 
equation, the value of the normal derivative of the function 
or of the function itself is specified along the boundary while 
the Cauchy problem for second-order partial differential 
equations involves the knowledge of both the function and a 
derivative. In the first two aerodynamic problems listed 
above, no absolute values are given but rather the jump in 
the value of the function along the boundary is prescribed. 
It is this type of problem with which the present report is 
particularly concerned. 

BOUNDARY- VALUE PROBLEMS IN PURELY SUBSONIC FLOW 

Since the purpose of this report is to extend the concepts 
of thin-airfoil theory used in subsonic theory to problems 
arising in supersonic theory, some discussion of the former 
will be given to provide lucidity as well as to furnish a basis 
for the analogy which will be shown to exist between the 
methods of solution arising in the two regimes of flow. 

The method whereby the solutions of the given problems 
can be effected is provided by Green’s theorem which relates 
a volume integral over a region R to a surface integral over 
the surface S enclosing R. If a, ft are any two functions 
which, together with their first and second derivatives, are 
finite and single-valued throughout R, then for the subsonic 
case 

5i> s — J / J W’“-W 8 < 7 > 

where the Laplacian operator, 

W2 = A 2 I 52 I — 
v dx 2 'dy 2 'dz 2 

is introduced and the directional derivatives on the left side 
are taken along the normal n , drawn inward, to the surface S. 
Identifying now the function a with the fundamental solution 
1/r and specifying that ft satisfies Laplace’s equation, equa- 
tion (7) simplifies to give 


( 8 ) 
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where 

i=[(x-x,) 2 + (y- yi ) 2 + (z-z, )*f 1 

The variables of integration in the equation are x h y i} z u 
while x , y, z are the coordinates of a point P either inside or 
outside of the region of integration. 

If the point P is assumed to be inside the region of integra- 
tion, it is evident that the function 1/r becomes infinite at P, 
and it is necessary to exclude this point from the region if 
formula (8) is to apply. Describing a spherical surface 2 
with radius e about the point P, and considering the integral 
over the two surfaces 2 and S which enclose the region, it 
can be shown that in the limit as e — >0 equation (8) becomes 


fixed and integrating over the low^er side of the xy plane, it 
follows that 


0=- 


4tt 




where the negative direction of the normal may be ignored 
since the integral is equal to zero. Subtracting these two 
equations gives the expression 


V(*,v, *) - - 4 ', | J,[0).(sf - a|)- 


( 10 ) 




(9) 


The physical significance of this last relation follows im- 
mediately: the term 1/r represents a fluid source and the 

term - represents a doublet with its axis lying along 

the normal to S, both source and doublet being situated at 
the surface point Xi, yi, Z\. The value of the function 12 at 
the point x , y , 2 is therefore given as an integral of a source 
and doublet distribution, the strengths of the two being 
determined directly from the respective boundary values of 12 
and 512 /dn. 

Equation (9) expresses the value of 12 in terms of the 
surface values of 12 and 512 /dn but this relation does not imply 
that a knowledge of both these variables is necessary for the 
determination of 12. As can be shown easily, another con- 
dition may be established which relates the two surface 
values. Applying equation (8) to the case where P lies 
outside the region of integration, it follows that the integral 
is equal to zero and that 12 and 5l2/5n on the surface are 
functionally dependent. 

Sufficient information is now at hand to provide a solution 
for the thin-airfoil boundary-value problems. Consider the 
region R bounded by the xy plane and a hemispherical dome 
of infinite radius lying above this plane. For all problems to 
which the results will be applied, the value of 12 may be as- 
sumed equal to zero at infinity. 1 The contribution of the 
surface integral over the hemisphere is thus zero and, from 
equation (9), 



the integral extending now only over the area r for which the 
integrand does not vanish. Equation (10) is the basic equa- 
tion from which all solutions in subsonic wing theory will be 
developed. It should be pointed out that the derivation 
proceeded from the assumption that the point (x, y, z) lay 
above the xy plane. When ( x , y, z) lies below the xy plane, 
however, the derivation can be carried through in exactly 
the same manner. Such a development reveals that equa- 
tion (10) is general and that no restriction need be imposed 
on the position of ( x , y , z) relative to the reference plane. 

As a particular application of equation (10) consider a thin 
symmetrical airfoil at zero angle of attack and set 12 =$ 
where <t> is the perturbation velocity potential. Conditions 
of symmetry demand that 

while 


and 


Thus, if w u —W[= Aw 

*=~lj jr AW 7 s dX ' dy ' (U) 

and the velocity potential is given by a distribution of source 
potentials. This distribution can be immediately related to 
the slope of the basic section by means of the equation 

w u _ _Wi == f dz\ 

V~ V \dx) u 


n u =4> u = ( S ) i=Ui 
512 W 


dx 

512, 

dz 


=w u 




where the integration extends over the entire plane and the 
subscript s indicates the function is to be evaluated at z x =0. 
The directional derivatives are necessarily in the direction of 
the positive 2 axis and subscripts are again introduced to de- 
note conditions on the upper side of the plane. Keeping P 


i If il is the perturbation velocity potential <l>. it is sufficient to assume that <t> and d$/c )n 
are zero on the sphere at all points having radius vectors which make finite (nonzero) angles 

with the positive x axis while $ and dfy'dx are merely bounded at all points infinitely distant 
from the lifting surface at a finite distance from the positive x axis. Thus, the results of the 
analysis can be applied to lifting-surface problems with their trailing vortex sheets. 


The symmetric airfoil can also be treated by replacing 12 
by the perturbation velocity w and in the case of the thin 
lifting surface with given loading the function 12 can be set 
equal to u. Employing, respectively, conditions of sym- 
metry and irrotationality, it follows that 512 u /dz — 512,/5z 
vanishes and, after setting A12 = 12 M — 12 , equation (10) becomes 

«(*’ V' Z)= L / Jr AQ (di, \) dX ' dy ' (12) 
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BOUNDARY-VALUE PROBLEMS IN PURELY SUPERSONIC FLOW 

Applications of Green’s theorem. — The problem to be dis- 
cussed at this point is the extent to which an analogue to 
equation (10) can be developed for supersonic flow fields. 
The first step in the presentation is, once more, the introduc- 
tion of Green’s theorem for equation (1) after it has been 
modified to the form given by equation (4). Employing the 
operator 

n 2 ==— 

dz 2 by 2 dz 2 

Green’s theorem now becomes 

- JJ (' lu) <»-/ X J rawcroffi (13) 

where u is the so-called conormal to the surface S and has 
direction cosines equal to v lf v 2 , v 3 such that 

vi=—n u v 2 =n 2 , v 3 =n 3 (14) 

where n 1} n 2 , n 3 are the direction cosines of the normal to the 
surface S (fig. 1). (The conormal at any point x h y 1} z x of a 
surface is the mirror image in the plane x=x x of the normal 
through the same point.) If <r and 12 are perfectly arbitrary 
functions, aside from satisfying the usual conditions of 
single-valuedness, etc., equation (13) represents an identity 
and this fact will be useful at a later time. For immediate 
purposes, however, <r and 12 will be chosen as solutions of the 
differential equation under consideration so that 

□V=D 2 12 = 0 

and, consequently, 

JJ ( 0 5;-'“) <iS - 0 < 15 > 

The use of equations (13) and (15) depends upon an under- 
standing of the physical nature of supersonic flow fields. 
The essential feature of such flow is the presence of Mach 
cones which correspond to the characteristic cones arising 
in the mathematical study of the wave equation. In ac- 
cordance with these concepts a disturbance in the flow field 
can affect the flow only within its aftercone, that is, the cone 
with vertex at the point of disturbance and with axis ex- 
tending in the direction of the undisturbed stream velocity 
vector; conversely, a point in the flow field can be affected 
only by disturbances which emanate from points within its 
forecone. 

When the disturbances are generated by a wing it is, 
moreover, necessary to speak more specifically about the 
nature of the leading edge of the wing. For all cases con- 
sidered here the assumption will be made that the plan form 
is a polygon, that is, is composed of straight line segments. 

If the wing is swept ahead of the foremost Mach cone, the 
cones arising at the leading edge will have as envelope a 
wedge-shaped surface passing through and extending back 
from the leading edge, while if the wing is swept back of the 
foremost Mach cone this cone will be the surface along which 
the air first experiences perturbations or disturbances. 
Thus, a point P with coordinates x , y y z is affected by all 


disturbances lying within its forecone T and at the same time 
behind the forward surface X, the nature of the latter surface 
being dictated by the leading edge. In figures 2 (a) and 
2 (b) these surfaces, along with the disturbance plane r, 
are indicated for two different wing plan forms. In the 
applications of equation (13) the volume integral is limited 
to the portion of space common to the surfaces T, X, and r 
and the surface integral involves a discussion of conditions 
on these surfaces. 

Up to this point the analogy between the subsonic and 
supersonic cases, insofar as the use of Green’s theorem is 
concerned, is quite apparent. The principal difference which 
occurred was brought about by the use of the true normal in 
the subsonic field together with the fact that the xy plane 
was covered by a hemispherical dome of infinite radius; 
whereas, in the supersonic field, the concept of the conormal 
was introduced and the volume to be considered was that 
enclosed within a finite region. In continuing the analogy, 
however, far more formidable obstacles arise. To begin 
with, the discussion of a and 12 over the surface in the sub- 
sonic case was relatively simple. Thus, with the exception 
of footnote 1, 12 could be assumed zero at infinity and A12 
was specified completely in the xy plane. But in the super- 
sonic case, although A12 can be assumed known in the xy 
plane and, as will be seen later, 12 may be evaluated on the 
forward boundary of the region, nothing is known of 12 on 
the forecone T. Hence a must be chosen properly so that 
the knowledge of 12 is unnecessary on T. The most obvious 
choice of <r would be a particular solution of equation (4) 
which would make <7=0 on T and this is in fact the choice 
used by Volterra (reference 1) and applied to aerodynamic 
problems in reference 2. However, if the analogy is to be 
maintained the choice of <7 is not arbitrary but must be the 
three-dimensional supersonic source corresponding to the 

z 



Vi s cos ( X-a ) 
v 2 ~ cos b 
Vj = cos c 

Figure 1— Geometric relations between direction cosines (m, n 2 , n j) of normal and direction 
cosines (ui, m, w) of conormal to surface S. 


Til r cob a 
n i -coB b 
n 3 = cos c 
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will tend to infinity as the surface S approaches r so that 
l/r c is not a regular solution to D 2 12=0 on T. However, 
as has been mentioned, equation (13) still must hold whether 
or not <t or 12 satisfy the wave equation and thus it still 
provides an equality. Hence, if the surface integral tends 
to infinity so also must the volume integral. Further, equa- 
tion (13) implies that these infinite portions just cancel since 
the difference of the two integrations must always give zero. 
To deal with such a problem quantitatively by the usual 
mathematical techniques would require the study of a limit- 
ing process for each new boundary-value problem. Hada- 
mard’s contribution was the introduction and justification 
of a concept which removed the necessity for studying the 
infinite portions involved. This concept is best presented 
by means of a new notation, thus the sign r is used and is 
to be read “the finite part of.” 

Using this concept it is possible to show that if a were set 
equal to l/r c , then equation (13) could be written 



(a) Triangular plan form. 

(b) Rectangular plan form. 

Figure 2.— Mach forecone from point P Or, y, z) intersecting surface t. 


fundamental solution 1/r in subsonic theory. But such a 
solution 

* =[(x— X,) 2 — (y— ViY— (2— Zi) 2 ] 

• c 

becomes infinite along the forecone r which has the equation 
(x—xi) 2 — (y—y i) 2 — (2— 2 i) 2 =0 

It is just this difficulty which apparently invalidates any 
furtherance of the analogy and the prediction in advance 
of an aerodynamic shape from a distribution of sources and 
doublets in supersonic flow. However, it is also precisely 
this difficulty which is overcome by Hadamard’s general 
methods. 

Extension of analogy by Hadamard’s method. —The full 
development of Hadamard’s methods cannot be given here, 
but a rough sketch of his reasoning will perhaps be useful. 
The basis of his arguments stems from equation (13). 
First it is admitted that the right-hand side of equation (13) 


inranf]- 


so that the “finite parts” of each side of the equation would 
be equal. Such a notation would, of course, in general be 
meaningless since in discarding arbitrarily a part which 
tended to infinity it would be possible, by proper combina- 
tions, to obtain as a remainder any finite value. The fact 
is, however, that the integrals involved in equation (16) tend 
to infinity only at a limit of the integration and this limit 
always involves the forecone T. It was consequently pos- 
sible to devise a manipulative technique to handle equation 
(16) without considering the singularities individually. It 
might be mentioned, without stressing the correspondence, 
that a treatment of improper integrals is also employed in the 
use of Cauchy’s principal value. In subsonic thin-airfoil 
theory and lifting-line theory integrals of the latter type are 
well known in the form 

MMdx o 
Jo X— X 0 


/i certainly tends to infinity as x 0 approaches x but the use 
of Cauchy’s principal value allows the very large values of 
the integrand obtained when x Q is on either side of x to just 
cancel in such a way that h is finite and unique. So again 
the integral 

T I f xc A(x)dx 

2 ~Ja (X 0 -X)*'* 

is finite and unique and given by Hadamard in the form 


or 


f '«A(z)-A(x 0 ) ^ 2A(x 0 ) 

a (z 0 -z ) 3/2 (Xo-a ) 1/2 


h= 



A' (x)dx 
tJxq — X 


2 A (a) 
■yjx Q — a 


It is actually possible to generalize the idea of “finite part” 
to the case when the exponent in the denominator is of the 
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form ./+1/2 where j is a positive integer but such a general- 
ization is not needed for aerodynamic applications and will 
therefore be omitted. 

In actual calculation, the evaluation of the integral I 2 can 
often be shortened. Thus, if the indefinite integral of 


A(x)dx 
(x 0 —x) 3/2 

is written in the form F(x)-\~C then 


/ 


2= lira T P 

X >xoLJ a 


A(xi)—A(x 0 ) 


dx i - 


2A(x 0 ) 


(x 0 —Xi Y n u '*' 1 (x 0 —a) 1/2 

It follows that if C is chosen so that 


< 7 = lira 

,X >J0 



then the expression for I 2 may be written I 2 = — \F(a.) + C] 
When C is chosen in this manner, the notation for the cal- 
culation may thus be modified to the form 


I A(x)dx~ r A(x)dx 
\Ja (X 0 -X )*' 2 Ja (X 0 — X ) 3/2 


and above the foremost disturbance. Applying equation 
(16) to the regions above and below the disturbance surface 
r (plane of the airfoil) yields the two equations 


f f rq 

d ( 

1 \ 1 df2~| 

J Jx+C+t[_ 

dv\ 

J‘c) r c du J 



-ffl" 


d 


dv' 

\r c J r c bv' 


dS=- 


\dS= 0 


(17) 


(18) 


where the prime indicates the surface value of ft on the 
opposite side of r from P. 

The integration over 2 can be computed for e very small. 
For convenience, consider P to be the origin; then it follows 
that since the conormal is in the X\ direction and the area 
element can be written 

y dy dd 

where 


and 


0=arc tan 


Vi 


y= vV+zi 2 


whore the asterisk indicates that the upper limit is not 
substituted into the indefinite integral F (x) + C. 

The technique for the calculation of the finite part has 
therefore been reduced to three simple steps: First, the 
indefinite integral F(x) + C is determined; second, the con- 
stant C in the indefinite integral is evaluated by means of a 
limiting process; and third, the lower limit of the integral is 
substituted into the indefinite integral and a minus sign 
prefixed. As an example, consider the integral 

l‘ x ° xdx xdx 

Ja (x 0 2 —X 2 ) 3/2 ~J„ (X 0 2 — X 2 ) 3/2 

In this case 

F ^ + C= (x 0 2 -x 2 y /2+C 

and 

r = h ™ [(2x 0 ) l ' 2 (x 0 ~x) l < 2 ~ (x 0 2 ~x 2 ) w*] =0 

so that, finally, 

X W-^y /2= ~ [^(a) +<71=^^172 

With the aid of this artifice the analogy between the sub- 
sonic and supersonic cases can be continued with relative 
ease. Thus, in equation (16) the right-hand member is zero 
provided we exclude the point P from the volume of inte- 
gration. This can be done most easily by limiting the 
integration to the x Y = constant plane, a distance e upstream 
from P. The portion of this plane intersected by the cone, 
and thus the section over which the integration must be 
carried, will be denoted by 2 (fig. 3). 

As drawn, figure 3 shows a cross section in a y x = constant 
plane for the special case when P is located directly behind 

843316—49 2 


the right side of equation (17) will give 


—lira T Q(x, y, z) ( dd | 6 ~ tydy 
c— »o LI Jo Jo (e — 


,2\ 3/2 "T” 


7 2 ) 


WEM f 2w dd 

fai Jo Jo V< 2 — 7 2 -l 


—lira 2irQ(x, y, z) ( , 6 ' 

«-»o L Jo (e — 

= 2tt^I(x. y, z) 


— eydy 


,2\ 3/2 


btt(x,y,x) f e ydy ~j 
dxi Jo Ve 2 -7 2 | 


7 2 ) 



Figure 3 —Cross section through region of integration used to obtain equation (19). 
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Hence the value of f2 at the point P, V(x,y,z), can be deter- 
mined from equation (17) with the restriction implied by 
equation (18). Further, since only the “finite part” is con- 
sidered, the integration over V yields zero and the two 
equations combine to give 

r,y,&= — 


When 12 is equal to one of the perturbation-velocity com- 
ponents, it is obvious that boundary conditions over X and r 
cannot be considered to be absolutely arbitrary since it is 
necessary to include the added restriction that the resultant 
potential 4> also satisfies the equation 

d 2 <k d 2 <t> d 2 <k _ 
dx 2 dy 2 dz 2 ~~ 



The only remaining difference between the subsonic solution 
for the distribution of sources and doublets, equation (10), 
and the supersonic solution, equation (19), is the integration 
over surface X. 


Considering first the case where the wing is swept behind the 
Mach cone, it follows that 


u(x i, y , z)dxi= ( v(x, 2/i, z)dy\ 
J Vz 2 +i' 2 J -V-r 2 -2 2 

= I w(x,y,zi)dzi 
J \x 2 -V 2 


and, after evaluating the partial derivatives of <$> and sub- 
stituting in the given differential equation, direct calculation 
leads to the conclusion that on X the following differential 
equations hold 


2 y 


dui 

w 


+2z 


dui 

dz 


+ Ui — 0 


DISCUSSION OF CONDITIONS ON SURFACE X 

By definition X is the surface on which the streamlines of 
the flow first experience pressure disturbances, that is, the 
surface along which the stream first becomes aware of the 
existence of the wing. Figures 2 (a) and 2 (b) were intro- 
duced to show the nature of the configurations involved for 
two different plan forms. It is apparent that when the wing 
is swept ahead of the foremost Mach cone the wedge-like 
form of X is comparable to the wedge appearing in purely 
two-dimensional problems while the wing swept back of the 
Mach cone has for its X-surface a cone and thus may be 
thought of as involving a purely three-dimensional problem. 

In order to determine the value of on X it is sufficient to 
impose the condition that the tangential component ol 
velocity is continuous across X. Such a condition represents 
no essential restriction since it is an immediate consequence 
of continuity of mass flow and continuity of the tangential 
component of momentum across the surface. As a result 
of this condition, however, it follows that the tangential 
component of the perturbation velocity is zero on the down- 
stream surface of X since it is obviously zero on the upstream 
surface. Moreover, velocity being equal to the gradient of 
the velocity potential the perturbation-velocity potential 
must be equal to a constant on X. But an arbitrary con- 
stant can be added or subtracted from the velocity potential 
so that with no loss of generality the value of on X can be 
assumed equal to zero and, since the conormal lies on the 
surface X, d<P/dv is also zero. 

The complete analogue to equation (10) has now been 
developed for &=$> so that 

tt(x,y,z) = 




0 dwi , 0 du\ 

2x bx +2y dy + Wl= ° 


where v h and Wi are the values of u, v, and w on X. The 
general solutions of these linear partial differential equations 
can be written as follows 




1 

' \lx 2 +y 2 



It has been stated, however, that the tangential component 
of the total perturbation-velocity vector vanishes on X, or, in 
analytical terms 

lui+mvi + nwi^O 

where l, m, n are direction numbers of any tangent to X and 
therefore satisfy the relation 

lx—rny—nz = 0 

Substituting the known expressions for u u v u it follows 
that 


l _ V 
y - \y 2 +z 2 





= 0 


or, using a different notation, 




Fo 




= 0 


Consider now the special case when Z=0 and rn= — nz/y 
Under these conditions 



(20) 




y 71 F 
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SO that 




Since the variables z/x and xjy are separated, the solution of 
this equation may be written 




K- 

y 


where K is a constant. Returning now to the case where 
Z, m, and n are in the ratio x: y : z , direct substitution into 
the above equation gives 




o 


so that 

- Fx 


y 

and 

f( 


+K^+-K -= 0 




This equation can, however, be written in the form 
F i(y/g) 


z/y+y/z 


-K- 


from which it follows that K= 0 and Fi~F 2 =F^= 0. All 
perturbation velocity components are seen to be zero, conse- 
quently equation (20) is valid for all cases in which the wing 
is swept back of the Mach cone. 

A discussion of conditions on the surface X will next be 
given for the case where the leading edge of the wing lies 
along the y axis (fig. 2 (a)) and where ft represents u or w. 
The perturbation-velocity potential $ may be given by the 
relation 

<t>{x,y,z)=^u{x u y,z)dx l 

where the plus and minus signs in the limit apply, respective- 
ly, to conditions above and below the xy plane. Since <t> 
must satisfy the basic differential equation, an added re- 
striction is imposed on u and as a result of this condition it 
can be shown that 

^4=0 

dz 


where u Y is the value of u on either the lower or upper surface 
of the wedge. It follows that the values of u on the two 
surfaces are 

Ui=fi(x,y) and Ui=^j 2 {x t y) 

and, since the solution is independent of 0 , and x is propor- 
tional to z in both cases, the final expressions are 

Ui=My) and ut=f 2 (y) 

If <f>(:r, y, z) is defined as an integral involving w, the same 
type of analysis leads to the conclusion that w on the two 
surfaces can also be expressed as functions of y alone. 
Perturbation velocity v will not be considered for this type of 


leading edge since the inclusion of u and w covers all com- 
monly used boundary conditions. 

It remains to substitute the results just obtained into 
equation (19) in order to study the contribution of the inte- 
grals over X. Apparently only one term in each integrand 
need be considered since the conormal is perpendicular to the 
y axis and the gradient of ft in that direction vanishes. As a 
preliminary step to setting up the integrals it is convenient 
to introduce a new coordinate system a;", y" , z " which is 
obtained by rotating the axial system about the y axis so 
that the x n and z n axes lie respectively in the lower and 
upper wedge planes wdiile the y" axis coincides with the 
y axis. The transformation of variables is 

(Xx — Zr) 

■yJZ 

z " = ~o ( x ' + 2 l) 

When tt=u, the last two integrals in equation (19) may 
now be written 


1 

2 7T 

1 I Cv+ V* 2 -* 2 

2tT I J y— ^x 2 —Z 2 


I fn-V * 2 -* 2 
Jy- V* 2 -2 2 


f c< 

fiiyjdyij^ 


f c< 

Myi)dyi\ 



dz"- 


dx" 


Substituting for r C) this expression becomes 


-a 

J.OI" 




y- Vr 2_ 2 2 


Myi)dyi 


dz' 


or 


V (*-72 I *")” 

f * d dx" 

Jo eta" 17 1 \2 7 r V 

VV _ V2 * ,/ )- ( *'-^ -(*+72^ ) 

_ 1 MVt)dyx 

2irjy-yg^g V x 2 —(y—yx) 2 —z 2 

l p+yg= ? jliyMyx 
2tt ]x 2 —(y—yxY—z 2 


It is apparent that if fi (y) = — jiiy) the integrals combine to 
give zero so that equation (20) may again be used in all cal- 
culations; moreover, the same condition applies when Q=w. 
The assumption that . f l (y) = — f 2 (y) is equivalent to postu- 
lating that in all cases ji{y) and j 2 (y) are odd functions of y. 
In application, however, this property is always maintained. 

It remains finally to consider the case when the leading 
edge of a wing is swept ahead of the Mach cone and when ft 
is a perturbation-velocity component. As a means of 
avoiding unnecessary complication in treating the problem 
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it is possible to substitute first the transformation (refer- 
ence 5) 

£= — f+mx 
r}= —x+my 
f= Iz 


This equation agrees with results given by Puckett (reference 
6) and others. 

As another example consider the solution used by R. T. 
Jones in reference 5 for a nonlifting symmetrical wing. Set- 
ting 9 equal to w in equation (20) and using the fact that w 
and u are related by the expression 


where the leading edge lias the equation y=mx , 2 = 0 , and 
Direct calculation shows that the basic differential 
equation and the Mach cone are invariant under the trans- 
formation and that in the new oblique coordinate system the 
leading edge lies along the r? axis while the planes of the wedge 
become 

m=o 

Because of the invariant properties of the transformation, and 
the fact that the 2=0 plane is fixed, equation (19) is appli- 
cable directly to the boundary-value problem for the swept 
wing m the new coordinate system. The treatment of the 
integrals over X can therefore be developed algebraically in 
exactly the same manner that applied to the previous case, 
hence u and w are constant along the lines 

£=Ff =0 , 77 = const 

and, again, if conditions of skew symmetry are maintained 
above and below the 2 =f= 0 plane the integration over the 
surfaces X cancel. Thus equation (20) is seen to be valid for 
9 equal to perturbation-velocity potential or perturbation 
velocity for all types of straight leading-edge configurations. 
And this is the complete analogue of the subsonic theory. 

APPLICATIONS 


f wdz 
ox J 0 


it follows that 


(22) 

For a wedge swept behind the forward Mach cone and having 
as the equation of its leading edge the relation y\=mx\ ) the 
expression for u may be written in the form 


Vd Z r dv A 

7 T Ox J VI 2_y2 Jo J Jy 


zdXi 


v_i [(x—Xi) 2 — (y~yi) 2 —z 2 f /2 


vliere 


Y 5-y-V( X “m) 2+22 (^ -1 ) 
* 1— 1 


77T 


5“1 


Performing the integration with respect to x if it follows that 

-( x ~m ) dy ' 


u= 


«s d p zdx r 

T bxJjjiZj, Jo 


[(y-yi?+^(x- y £) - (y-yiY—i 


INTERPRETATION OF PREVIOUS RESULTS 

As a means of indicating the various problems to which 
equation (20) can ho applied, three previously published re- 
sults will be discussed. These applications include, first, the 
expression for the perturbation-velocity potential of a sym- 
metrical nonlifting airfoil (reference 6), second, the calcula- 
tion of pressure distribution over a semiinfinite wedge with 
leading edge swept back of the foremost Mach cone (reference 
5) and, third, the integral equation method for determining 
the load distribution over a lifting surface of arbitrary slope 
(reference 7). 

As in the case of equation (11) for subsonic flow, let 0 
represent velocity potential < i> and consider the case of a sym- 
metrical wing at zero angle of attack. 

Then dft u /c) 2 =w K and b9.,/dz=w h where w u and w, are in- 
duced vertical velocities on the upper and lower surfaces, 
respectively. Moreover, <i>,=0 for the symmetrical 
case so that, since w u — Wi=2w u , 



7T 



dxidyi 


The integral in this equation is finite at T so the finite part 
sign may be disregarded and 


<f>= — 



w 


- dxrfiji 


and, after reversing the order of integration, 
_w u d (' Yi , f 2 

"~ T ado 

-{ x ~ V m) zdz 


-In 


[(2/ — 2/i) 2 +^ 2 ]-y/ ( y—y\Y z 2 

=C? P- 

27ro.rJ„ 


_V( x ~S) 2 “ (2/_yi)2_22+ *( a:_ S)_ 


dy 1 


where k= 1 for 0 ">)/ 1 /m and —1 for a:<y,/m. Taking the 
partial derivative with respect to x and noting that the value 
of the logarithm at the upper limit is zero, the value of the 
induced velocity is 

= _w u t‘ Y ' _ _ dy, 

and integration yields the final result 

m , (a— my) — V (1 — m 2 ) (x 2 —y 2 —z 2 ) 

" ir A /l— m 2 ^(mx—y) 2 -\-z 2 (l—m 2 ) 


(21) 
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Using equation (2) and setting w u /V equal to (< dz/dx) 0 the 
slope of the surface, this may be written 

q = — 2 (dz\ m ln (x— my} — V(l — m 2 )^— y 2 — z 2 ) 

” * \dxj o ^Jl—m 2 V {mx—y) 2 -\-z 2 {\ — m 2 ) 

(23) 


Equation (23) gives the pressure coefficient at any point in 
the field produced by a wedge swept behind the Mach cone. 
When 2 is set equal to zero the pressure distribution over 
the wedge itself is determined and the equation corresponds 
exactly with equation (12) of reference 5. 

When loading is to be prescribed over a thin lifting surface, 
Q may be assumed equal to the perturbation velocity u. A 
direct consequence of this assumption is that in equation (20) 

bQ 

bz ~~ dz 


since, from conditions of irrotationality, 


By definition 


bu u _bw u _ bw t _ dui 
bz bx ~~ bx ~~ bz 

7 2 Uy 2l/ l 

y ’ ° p — y 


and load distribution in coefficient form A p/q is given by 
the relation 

A p Cv u 

so that 


Ap 

T 


2ui . 2u u _ n Au 
y^~ y~ z y 


Equation (20) can therefore be written in the form 



= J I f Xp zdx l dy l 

" 4*1 JJ r q Vt (x- Xl ) 2 -(y-yij 2 -z 2 ] 3 

If equation (24) is transformed to the original 
variables, the relation for u is 


(24) 

space 


u(X, Y, Z) 


Ap 


ZdX x dY l 


P f f g 1 1 


) 2 +Z 2 ]} 


3/2 


(25) 


Equation (25) is valid for arbitrary plan forms with known 
load distributions. Particular examples which may be 
worked out with relative ease are the lifting surfaces carrying 
constant load. Once u is known the value of w can be 
determined from the integral 


b f* 
W dz)-* 


udx 


and from w the ordinate z of the surface as a function of x 
and y is given by 

w 




v 


dx o 


where l. e. denotes the leading edge. A discussion of 
trapezoidal, rectangular, and triangular plan forms with 
constant loading is given in reference 2 although the method 
of derivation is different. 

Interest in constantly loaded wings has been based pri- 
marily on the fact that in certain cases they can be combined 
to produce surfaces of given camber. Thus, a superposition 
of trapezoidal plan forms of variable rake, the constant 
loading over each trapezoid being a function of its rake 
angle, can be used to produce a flat plate of trapezoidal or 
rectangular plan form at an arbitrary angle of attack. In 
this case the loading as a function of rake angle is determined 
so that induced vertical velocity is kept constant. For 
problems in conical flow a lifting element can be constructed 
by subtracting from a constantly loaded right triangle with 
angle of sweep equal to 8 the constantly loaded right triangle 
with sweep angle equal to 8-d8. The resultant element carries 
a constant load and has a sweep angle equal to 8. By sum- 
ming these elements it is possible to find the load distribution 
as a function of 8 such that certain flat lifting surfaces at 
angles of attack are formed. In reference 2, triangular 
wings swept back of the Mach cone were studied by this 
method for arbitrary angles of yaw. Brown (reference 7) 
has used this same lifting element to study the more restrictive 
case of the symmetrical triangular wing. 

A brief discussion of differences existing between the meth- 
ods for producing the swept-back lifting element will shed 
some light on the various lines of attack. The approach 
used in reference 2 is essentially mathematical in that a par- 
ticular solution of the partial differential equation is used in 
conjunction with Green's theorem to satisfy the boundary 
conditions of the problem. The principal criticism of such 
a method is that the physical interpretation is missing. The 
use of equation (25), however, removes all such criticisms 
for precisely as in the case of incompressible flow the lifting 
element is created by distributing doublets over the wing. 
In Brown's solution it was necessary for him to determine 
first a line of sources by means of an integration along the 
line and then to form the doublet line by differentiating along 
the normal to the line. The order of differentiation for in- 
compressible flow is immaterial, since the limits in the inte- 
gral are independent of the position of the point P at X, Y, Z. 
Supersonic flow destroys this property and it is only after the 
introduction of the concept of “finite part" that the deriva- 
tive of an integral may be written as the integral of the dif- 
ferential coefficient of the integrand. Equation (25) thus 
maintains the analogy with previous work. 

LOAD DISTRIBUTION FOR ROLLING WING 

The usefulness of equation (20) is not at all restricted to a 
synthesis of previously known solutions. As an example of 
its generality consider its application to the problem of the 
rolling wing with leading edge swept ahead of the Mach cone. 
Figure 4 shows the boundary conditions involved. The value 
of w is specified over the wing and, since the Mach cone is 
behind the leading edge, the value of the perturbation veloci- 
ties u, v , and w are of course zero ahead of the leading edge. 
Assume for the moment that a symmetrical body at zero 
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angle of attack is considered. It follows that if tt v —w u and 
Q t = Wi then equation ( 20 ) can be written in the form 




(26) 


since, for reasons of symmetry, the normal gradients of w on 
the two surfaces are equal. Using now the fact that the 
Mach cone is behind the leading edge then the pressure over 
the upper surface is independent of the shape of the lower 
surface and equation (26) may be applied directly to the roll- 
ing flat plate if w u is determined from the given induced 
velocity on either the upper or lower surface. This method 
of approach, of course, limits the solution to cases where the 
leading edge is ahead of the Mach cone. 

If the rate of roll is given as P radians per second then 
2w u = < 2PY\ and equation (26) becomes 


w 


P|ff P 2 YiZdYidXi 

7T J JrKX-xy-piY-Yty-pZT 2 


(27) 


The area r in equation (27) is that contained between the 
leading edge and the trace of the forecone on the XY plane. 
Figure 5 (a) shows the configuration for three traces corre- 
sponding to forecones from the points Pi, P 2 , and P 3 . The 
region containing the point P 2 is distinguished from that 
containing Pi and P 3 by the fact that r for P 2 lies ahead of the 
Mach cone from the apex and, furthermore, entirely on the 
right of the X\ axis. The regions corresponding to Pi and 
P 3 differ in the fact that when integrating from + <» to Z 
to find u, the upper limit of the integral in the first case is 
the Mach cone X 2 — fPY 2 — /5 2 Z 2 =0 whereas in the latter case 

mX—Y 

the upper limit is the leading-edge wedge (fig. 

5 (b)). The solution must be carried out separately for each 
of these regions but only the details for the region corre- 
sponding to Pi need be given here since the others are similar. 

It follows that the induced velocity u at the point P is given 
as the sum of the triple integrals 


5 Pp f° V j V C* &Z4X, 

~5XxJ. V * 5z ^/ Z °h(i.z» Y ' d) J_V, {(X-X^-pHY-Yy-fZ, 


!]3/2“ 


d Pfz /y r* /3%rfZ, 

dX 7r Ji V ^_ WJ dA 'J 0 1,d .]> n l [(X-X.r-dYY- Y^-pZoT 2 

_ s j>_rr* z,r,(x +t $>r, 

bXir Jl t / x >- py > ° JMk.z*) [(T-Y{) a +ZJ]J(x+k ^‘ ) — 18 2 ( F— Yi) 2 —p 2 Zo- 


(28) 


where A{k,Z 0 ) is the value of Y x determined by the intersection of the forecone with the leading edge on the right and 
left sides with k equal respectively to — 1 and + 1 . (See fig. 5 (a).) Thus 


A(k,Zo) 


-(* |+yr)+wV(-g+* iQ+a.’Q?-?) 


in 


-P 2 


After reversing the order of integration and integrating with respect to Z 0 , it follows that 

p c) r° 

S Y4YJn 

A = — 1,1 ^7T OA. J A (k,Z) 


' (*+* C-)-V(r+<:S)'-g , (r-r.) , -yg 

_( x + k S’) + V(‘ Y+t 


Moreover, since the integrand is zero at Y y =A(Jc,7 . ) the derivative with respect to X can be taken inside the integral and 

^ P 7 f 0 M 

I=-1, 1 7T Jaoc, Z) ^ x+ Y k j_ j32( Y^ Yl y- p-Z 2 

Integrating in this equation and combining terms it follows that induced velocity u is given by the expression 


P . mpY+X 


100*7+ X 


U ~1T m2 (m 2 ^ 2 -l) 3/2 i arc 81 " [uV(mZ+r) 2 +(l 


— m 2 0 2 )Z 2 j 2 ) 7 r 


m 


mf¥Y—X 


r 


mpY-X 


(m 2 /3 2 — l) 3/2 j aiC Sin U V(mX- F) 2 + (l-m 2 0 2 )Z 2 ] + 2 

(29) 


Setting Z =0 in equation (29), pressure coefficient C„= — 2 u/V is given by the expression 


„ 2 mB 2 Y—X . |~ mp i Y—X~\ 

Cv ~TrV m ~\ (m 2 /3 2 — 1) 3/2 arC Sm l_j3(mZ— Y) J 


mp 2 Y+X 
(m 2 \ 3 2 — 1) 3/2 


9 arc sm 


f m/5 2 Y + X~] . mfYYit 


[J(mX+ Y) J (m 2 /5 2 — 1 ) 


3/2 


(30) 


when A')>/3 
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Sect ion S-S 


Sec i/on A - A 



(a) Planform. 

(b) Sections showing distribution. 

Figure 4.— Vertical velocity distribution for rolling wing. 




Figure 6.— Variation of pressure coefficient for rolling triangular wing swept ahead of Mach 

cone. 


Z x 

| ,-Mach cone, X z -ft 2 Y 2 ’P z Z* O 



(a) Traces in Zi=0 plane. 

(b) Traces in Xi= constant plane. 

Figure 5.— Integration regions used in determining loading over triangular wing swept 
ahead of Mach cone. 

This solution holds in both regions containing the points 
Pi and P 3 . However, in the region ahead of the Mach cone 
but still on the wing (region corresponding to P 2 ) it is easy 
to show that 

2jP 771 ? 

C v = V(m 2 p 2 1) 3/2 (^-*) (31) 


where Ylm<CX<CpY. Figure 6 shows a spanwise plot of 
C P £l/ (J^) J for m= 2 and 0=1. 

Equations (30) and (31) provide sufficient information for 
the calculation of the stability derivative for damping in roll, 
C lp . Integration of the load distribution yields the result 
that 

p dC 1 1 

^“d(P6/ 2V)~ 30 

LOAD DISTRIBUTION FOR PITCHING WING 

Another simple application of equation (20) is found in 
the solution to the problem of the pitching wing. Figure 7 
shows the boundary condition involved which is that the 
vertical induced velocity be a linear function of X t . If Q 
is the rate of pitch in radians per second, then w=QX x on 
the wing. Again the solution is obtained only for wings 
which have leading edges swept ahead of the Mach cone. 
(Although solutions can be obtained for leading edges swept 
behind the Mach cone, they involve integral equations and 
do nothing to illustrate further the direct methods of this 
report.) 
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(b) 


(a) Planform. 

(b) Sections showing distribution. 

Figure 7.— Vertical velocity distribution for pitching wing. 



Figure 8.— Variation of pressure coefficient for pitching triangular wing swept ahead of Mach 

cone. 


Ill the rolling wing ease, £2 was set equal to perturbation 
velocity w and as a result a distribution of doublets was 
used in equations (27) and 28). As an example of the 
manner in which source-sink distributions may be used for 
the same type of problem, equation (22) will be applied in 
the present case. Since the wing is swept ahead of the fore- 
most Mach cone, induced effects on the upper and lower sur- 
face are independent and 


1 f f QX x dX x dY x 

ttJ Jr V(AT-Ar l ) 2 -/3 2 (F-r l ) 2 -^Z 2 


(32) 


Again three regions containing the points Pi, P 2 , and P 3 are 
distinguished (fig. 5 (a)) and the solutions will be derived 
only for the region containing Pi. Integrating first with 
respect to Y y and then differentiating with respect to X 
yields 


b<t> Q , v 

u== bX = - p {X ~ 

Y—kmXi 


■JiX-Xtf-pZ 2 


PZ) + ? B 2 

irp t = -i, ! 

dX x 


; -X 


B(k,Z) 


Xi ^ y arc sin 


(33) 


where 

_ 7 _ X— & 2 Ykm — fik V (Xm —kY) 2 -\- Z 2 ( 1 — ra 2 /3 2 ) 

B(k,Z)=m 


Considering the limit as Z— > 0 and integrating gives: 


C p X q (m 2 /3 2 —l) 3/2 =~ V (X 2 —/3 2 Y 2 ) (m 2 /3 2 — 1) — 


mX{2— m 2 f} 2 ) — 


mX{ 2 — m 2 /3 2 ) — Y 

7T 

mX(2-m 2 l3 2 ) + Y 

7T 


arc sin 


in( 


arc sin | 


X-p 2 Ym\ , 
pY—f3mXj' 

X±fYm\ 
/3F+ X ) 
(34) 


Formula (34) is valid for (he regions P x and P 3 of figure 5 (a). 
For the region P 2 the solution is: 




V 
2 Q 


(m 2 p 2 — 1) 3/2 = Y-2Xm + Xm 3 p 2 


(35) 


Figure 8 shows a spanwise plot of ^ j^ 1 / * or 
and /3= 1 . 


m- 


Equations (34) and (35) provide sufficient information for 
the calculation of the stability derivative for the damping in 
pitch Cm Q . Integration of the load distribution gives the 
result 

dC m _ 4 

d(Qc/2V) 9/3 


where the axis of rotation is at X~2cj?>. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., January 1948. 
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Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Designation 

Sym- 

bol 

Longitudinal 

X 

Y 

Z 

Lateral.... 

Normal 



Force 
(parallel 
to axis) 
symbol 


X 

Y 

Z 


Moment about axis 

Angle 

Velocities 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 

Rolling 

L 

Y >Z 

Roll 

6 

u 

T1 

Pitching 

M 

Z >X 

Pitch 

'r 

0 

v 

r 

a 

Yawing 

N 

X >Y 

Yaw 


tv 

r 


Absolute coefficients of moment 

fy /i A/ 

Vl ~qbS Cn ~^S 

(rolling) (pitching) 


C n = 


N 


'qbS 

(yawing) 


Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript. ) 


D Diameter 

p Geometric pitch 

pjD Pitch ratio 

V' Inflow velocity 

V, Slipstream velocity 

T Thrust, absolute coefficient C T - 


PROPELLER SYMBOLS 

P 

C, 


Q Torque, absolute coefficient C Q 


1 hp=76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower=0.9863 hp 
1 mph= 0.4470 mps 
1 mps=2.2369 mph 


'p7l 2 Z> 4 

Q 


v 

n 

$ 


~ prPD 6 

5. NUMERICAL RELATIONS 


X 

Power, absolute coefficient C P = — 

pn 3 D s 

Speed-power coefficient =-^p ^2 
Efficiency 

Revolutions per second, rps 
Effective helix angle =tan -1 ^~^) 


1 lb=0.4536 kg 
1 kg= 2.2046 lb 
1 mi= 1,609.35 m=5,280 ft 
1 m=3.2808 ft 


